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A3STBACT 

The  theory  of  the  three-dioensional   rotational  flow  of  an  In- 
conpreasible  and  inviseld  fluid  through  an  axial  tur"bomachine  is 
described  and  the  hydro dynanical  equations  are   simplified  by  consid- 
ering an  infinite  number  of  blades  in  er^ch  row.     The  forces  of  the 
blades  on  the  fluid  ere  treated  as  non-consenr.^tive  bocly  forces  dis- 
tributed uniformly  about  the  axis. 

FomiulFition  of  the  mathemp.tic?OL  problem  le^ds  to  one  non-linear 
partial  differential  equation  and  t/o  integral   equations  for  the 
three  Telocity  com;)onents,     A  linearized  solution  of  these  slaultan- 
eous  equ  tions  for  any  prescribed  bl.'sde  loading  is  based  on   the 
consider/: tion  that  the  vorticity  .^^ener  ted  by  the  blades  is  trans- 
ported do'jrnstpea.'n  by  the  momx  axial  Telocity,     An  iteration  process 
which  leads  to  solutions  of  greater  accuracy  ie  deTeloped  by  consid- 
ering for  each  iteration  that   the  Torticity  is  transported  by  the 
Telocities  found  ty  the  preTious  iter-tion. 

The  Bessel's  functions  w.ich  occur  in  the  Green's  function  solu- 
tion are  replaced  by  their  asyinptotic  Talues  and  the  infinite  series 
is  suniued  to  express  the  solution  in  closed  form.     The  iter -tion 
process  is  then  adapted  to  nechanical  cjJ.culations  by  diTiding  the 
region  of  Torticity  into   snrJ.i   rings  of  rectangular  cross-section  and 
detorTiiining  the  influence  on  the  velocity  of  a  unit  change  of  vortic- 
ity  in  e  ch  of  thene  riu'^s.     Once  this  influence  is  eatf^bliahed  it 
is  relatively  easy  to  calculate  the  velocities  in  any  axial  flow 
m-^chine  with  any  T;re scribed  bli^de  loading. 
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IHTHODUCTIOK 

The   derelopment  of  tarboaachinezy  until  recent  years  has  resulted 
primarily  from  empirical  nethods,  both  in  design  and  analysis,   Bscently 
theoretical  siethods  have  led  to  the  adranceaent  of  axial  flow  ooinpressors 
and  turbines  eren  thou^,  because  of  the  complicated  matheisatieal  problem, 
many  sistplifyin^  assunptions  vero  necessarily  made  and  solutions  were 
ap-^roxinate. 

In  analogy  vith  the  three  dimensional  ving  theory,  problems  of 
flow  in  turbomaehines  may  be  classified  as t 

(1)  The  Direct  FroblsBi  The  direct  problem  of  calculating  the  flow 
in  turbomachinexy  is  th^t  of  detemining  the  velocity  field* 
the  blade  forces,  and  the  distribution  of  energy  in  the  fluid 
when  the  blade  shape,  the  blade  speed,  end  the  appropriate 
boundaxy  conditions  are  prescribed. 

(2)  The  InTerae  Problem!  The  inverse  problem  of  calculating  the 
flov  in  turbomaohiaery  is  that  of  detenining  the  velocity 
field,  the  bl&le  shape,  and  the  distribution  of  enexo'  i&  ^h« 
fluid  vhen  the  blade  loading,  the  blade  speed,  and  the  boundary 
conditions  are  prescribed, 

(3)  The  Intermediate  ;'roblep:  An  intezvediate  problem  which  ap- 
pears to  be  of  interest  is  that  of  determining  the  velocity 
field,  the  m/ignitude  of  the  remaining  forces,  and   the  distrib- 
ution of  energy  in  the  fluid  when  one  blade  force  is  prescribed 
and  the  blade  shape  is  partially  prescribed. 
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The  direot  problen  arises  when  it  is  desired  to  inrestigato  a 
giren  machine  "off  the  design  point". 

The  inverse  probla-n  may  be  solyed  for  the  initial  design  but 
usually,  because  of  the  st2*uotural  limitation  on  blade  shape  (see 
appendix),  it  will  be  more  practical  to  formulate  and  solre  the  inter* 
mediate  problem* 

In  order  to  simplify  the  mathamatioal  problem  the  fluid  will  be 
assumed  to  be  invlsoid  and  inoompreseible  and  the  blade  forces  will 
be  treated  as  body  forces .  uniformly  distributed  throu^  the  fluid* 
so  that  the  flow  is  symaietrical  about  the  axis  of  rotation  and  the 
vorticity  is  no  longer  shed  in  sheets  behind  each  blade  but  is  contin* 
uously  distributed  orer  the  region  downstream  of  the  blade  row. 

An  axial  flow  machine  for  which  the  inner  and  outer  boundaries  con- 
sist of  concentric  circalar  cylinders  extending  to  infinity  in  the  di- 
rection of  the  flow  will  be  conaidered. 

The  niathematical  problein  is  formulated  by  considering  the  tiiae 
rate  of  channel  alon^'^  e  atreamsurfece,  of  the  tancential  component  of 
the  vorticity  vector.  The  diffic  Ity  of  this  problem  lies  in  the 
solution  of  the  non-linear  partial  differential  equations  that  describe 
rotational  fluid  motion.  In  order  to  overcome  this  difficulty  a  method 
of  iteration  is  developed  whereby  solutions  of  any  required  degree  of 
accuracy  may  be  obtained. 

The  first  step  of  this  iteration  j^^ocess  provides  a  linearized 
solution  (cf .  tiSarble  )  based  on  the  assumption*  analogous  to  the 
Prandtl  three  dimensional  wing  theory,  that  the  vorticity  is  •trans- 
ported" downstream  by  the  mean  axial  velocity  and  is  not  influenced 
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by  its  ovn  Induoed  Telooitlat*  This  linearized  aoluUon  is  not  as  ne- 
curate  as  ni^t  be  expected  from  the  abore  a&aloigr»  one  reason  being 
that  the  Yorticit7  is  shed  in  three  dimensional  spjice  instead  of  in 
a  two  dimensional  sheet  so  thnt  the  induced  velocities  are  likely  to 
be  reiy  large. 

?or  the  second  step  of  the  iteration  process  the  Torticity  is 
considered  transported  bgr  the  velocities  found  by  the  linearised  so- 
lution of  the  first  step.  The  succeeding  approximations  are  obtained 
in  the  same  imy,  in  each  case  using  the  velocities  of  the  preceding 
approximation. 

The  solution  is  obtained  bgr  finding  the  appropriate  Green's  func- 
tion -  th&t  is,  a  function  Q  (r,  z;*<»^  )  which  gives  the  velocity, 
consistent  with  the  boundaxy  conditions,  induced  at  any  point  of  a 
circle  r,  z  by  a  unit  change  in  tangential  vorticity  of  a  vortex  ring 
at  radius  o<  at  an  axial  coordinate  /3  •  The  Besoel  functions  which 
arise  in  this  solution  ars  replaced  by  their  asymptotic  values  rJid  the 
Green* s  function,  vhich  would  occur  as  an  infinite  sum  of  Bessel  func- 
tions, is  l^  this  means  expressed  in  closed  form. 

The  resulting  expression  for  the  radial  velocity,  a  double  inte* 
gral,  is  modified  and  adopted  for  mechanical  calculation. 
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I,  notation  and  Symbolc 
Itie  flow  is  dosoribed  (?if;>  1)  in  a  o/lindrioal  coordinate  aystea 


FttP-  t   2t  by  the  velocity  components 
u,  T,  w,  respectively.  The  corres- 
ponding radial »  tan^eatial,  and  axial 
velocity  components  are 


S    = 


0  z. 


5  w 
Or 


".I 


ne.  1 

Coordinate  CysteiH  and  Desig- 
nation of   /elocity  and  7ortio« 

ity  Co.T.poneats 


In  addition  we  will  use  the  following-  synbDls. 


cO  =  angular  velocity  of  rotor 

p  -  pressure 

V  -  absolute  velocity  vector 
n   -  vortioity  vector 

F  =  force  vector  ^ force  of  blades  on  fluid) 
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II.  ^jrdrodynainio  Equations 
The  folloving  equationB  are  simplified  far   this  case  of  an  invlsoid 
and  incompressible  fluid  in  stead/t   adiabatic,  azially  sycinotriCQl  flow* 

Sg^gations  of  llotiont 


Or  r  O^  r       Or 


r    dr  oi 


7> 


0^  -^i^'^./r 


Qoatinulty  liquation* 
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III.  7ozmuIation  of  the  Stothematloal  Problem 
a)  Developtaent  of  the  Basic  ZauatlonB 

iieoause  of  the  axial  syionetry  only  the  tangential  rortioity  is 
associated  with  the  radial  and  axial  relooitles*  while  the  radial  and 
axial  Torticity  ooraponents  are  only  asaociated  with  the  tangential 
velocity*  The  tangential  vortioity  constitutes  an  annular  Tortex  ring 
and  by  considering  the  defarmation  of  this  ring  information  concerning 

the  radial  and  axial  velocities  can  be  gained*  It  is  well  known  (of* 

2        1 
Ueyer   •  Marble  )  that  the  circulation  about  a  defondng  vortex  ring  is 

constant  in  a  conservative  force  field*  Only  a  radial  stretching  can 

occur  since  the  flow  is  aymsetrical  about  the  axis  so  that  the  constancy 

of  ciroalation  requires  that  in  a  conservative  field 

•  f?  (?)  -  °  w 

In  the  presence  of  non- conservative  blade  forces  or  if  radial 

vorticity  is  present  the  circulation,  and  hence  the  quantity  *^,  will 

r 

not  remain  constant*  The  law  roveming  its  variation  can  be  derived 
froQ  the  equations  of  motion*  The  equations  of  lootion  in  vector  fona 
can  be  written 

l7  X  i2  ^  -^  ^/^  ^  i^^^-  r    -  ^  (2) 

If  we  take  the  curl  of  both  sides  and  use  the  equation  of  continuityi 


cr  •  ^  =  o> 


(3) 
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thar«  results 

^o  tangsntial  conxponent  of  this  vector  aqiL-tion  it 

^  fr  f  *^ yr  -^  -;r     '^  ^r  ^  ^  91.   -^  -^   -f-  sir     ir?^ 
which,   CLsin^  the  definitions  of  Torticity,   can  ^e  sioiplified  to 

Sinoe  the  circuit  tion  is  directly  pro.jortional  to  the  quantity—  , 
this  equation  expresBes   the  Iav  governing  the  tine  ch n^e  of  circula* 
tion  around  the   tangential  rorticitj  ring.     The  first  terra  on  the  right 
is  the  oh<mge  in  the  axial  direction  of  the  centrifugal  force,    the  8eo-> 
ond  term  is  the  chazige  of  the  radial  blade  force  in  the  axial  direction, 
and  the  third  term  is  the  change  of  the  axial  blade  force  in  the  radial 
direction.     A  little  thou^t  will  show  that  these  terms,  in  each  case, 
represent  moments  tending  to  cause  rot -tion  of  a  particle  about  a  tan- 
gential axis.     It  is  in  thi3  mpjmcr  that  a  non-oonserratire  force  field 
tends  to  effect  n  chan,^©  in  the  tangential  rortex  ring. 

The  tine  derlTatire  on  the  left  is  taken  along  a  streaialine  and 
is  written 
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In  order  to  compare  the  relative  magnitude  of  the  terme  on  the  right 
we  uae  the  definition  of  tangential  velooity  and  write 

Hemembering  that  for  the  axial  flow  maohine  the  radial  velooity  is 
onall  oocapared  to  the  axial  Telocity  and  expeotin^  from  physical  oon* 
siderations  that  the  velooity  distribution  will  be  smooth,  it  appears 
that 

"nils  inequality  will  be  useful  for  the  first  approximations  but  must  be 
more  closely  investigated  for  the  final  api^roxiaatione.  It  is  conceiv- 
able that  for  flow  that  differs  greatly  fron  vortex  flow  the  radial 
ohnnso  of  taT^ential  vortioity  would  be  of  such  magnitude  that  this 
inequality  is  not  justified.  It  will  however  be  useful  to  group  the 
smaller  terrpjg  separately,  that  is,  to  write 

:^  -  ^ l:^^~r^  -^ ^ir     TT '"^JT^  rln     ^-fIJ  v7) 

Using  the  definition  of  vortioity  we  can  write  an  equation  for  the 
radial  velooity  as  follows 
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ao  that  the  left  sidot   if  set  equal  to  zero,  would  be  a  linear  partial 
differential  equation  for  vhioh  a  solution  is  Icnown. 

Since  the  fluid  is  inrisoid  the  force  of  the  blade  on  the  fluid 
will  act  nomal  to  the  blade  and  hence  noxroal  to  the  relatire  velocity 
ao  that 


u 


/;.    ^  f^-<^r]  f^     .-    ^  /^z     =   ^  (9) 


The  radial  velocity  and  the  radial  force  are  both  anall  compared  to  the 
other  force  and  velocity  components  so  that  with  very  good  approxima- 
tion 

^z   -  -  ^ P-^  (10) 

Using  this  relation  in  the  equations  for  radial  velocity  we  have  one  of 
the  final  equations  of  the  iteration  process 

For  detenninin^  a  relation  between  the  tancontial  velocity  and  the 
tansential  force  we  have  the  second  equation  of  laotion 

which  states  that  the  variation  along  a  streamline  of  the  raoojent  of 
moraentun  is  equal  to  the  mcnent  of  the  tancentiel  force  about  the  axis 
of  rotation.  Then  since  we  will  prescribe  the  tangential  force  we 
oan  express  the  tangential  velocity  as  an  integral  along  a  cylinder. 
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/3=  - 


The  continuity  equation,  when  integratedt  prorides  an  equation  for 
the  axial  velocity  in  terma  of  the  radial  yelooity 


^         —-         L--^    ^  ^^  (13) 


We  have  developed  three  baaio  equations,  Qqa*  11,  12,  13 t  to  be 
used  in  formulating  the  linearized  problem,  for  the  first  approximation, 
and  in  the  oonatruotion  of  an  iteration  prooedure  whereby  more  exact 
approximations  t^^   be  obtained* 
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^)  "Rie  Uxiearlzad  Problea 

In  order  to  lixieariza  tho  bafiio  e(iuatloiiB  derired  in  the  prerious 
aeotion  wa  shall  consider  axial  flowi  where  the  radial  velooity  is  very 
amall  oompared  to  the  axial  Telooity* 

u  <  <  w 

and  will  aaanne  that  the  rortioity  la  transmitted  axially  downatream 
with  the  mean  axial  velocity  i|  ao  that 


^.(-^}  =  ^^ii(^) 


furthermore^  if  the  radial  force  is  onalli  we  have  for  calculating  the 
radial  velocity  the  expreaaion 


=  ;;7/o-|/^^V^^^^^^W/  (H) 


For  calculating  the  linearized  tangential  velocity  Eci*  12  becooea 
and  the  linearized  axial  velocity  ia,  from  'Cfi»   13 1 

We  are  now  able  to  outline  a  solution  to  the  linearized  problem* 
We  will  consider  here  the  inverse  problem  and  will  specify  the  two 
force  components  1^   and  F^  .  ^  will  prescribe  Fj.^-^  ?^  and  F^^  l^(r,z), 
and  will  specify  the  bo'jndary  conditions  applicable  to  axial  flow. 
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First  tha  tangential  velooity  is  determined  by 

V    -  /'.^  -'  ' —  ^/^  (13) 

Knowing  the  tangential  velooit/  we  oan  determine  the  radial  relocity 
froi^  2(1«   14 


e>  r  »-  "'' 


with  the  boundar  '  condition  that 

Ve  have  here  an  olllptio  partial  differential  equation  for  tlie  radial 
Telocity.  Since  we  know  the  right  hand  sidOf  the  unhomogeneous  parti 
and  hare  complete  boundary  Talues  we  oan  solve  ttiia  equation  for  the 
radial  velocity.  ITiis  equation  was  solved  by  Marble  in  terms  of  a 
Green's  function  GCrtz;  ^,fi)   which  gives  the  velocity  induced  at  any 
point  of  a  circle  T%t   by  a  unit  chan^'e  in  tangential  vortioity  of  a 
vortex  ring  of  radiua  o^  at  an  axial  coordinnte  /3  .  By  Marble's  solu- 
tion the  radial  velocity  is  exxjressed  as 


u. 


where 


^     J         J         ^^-<.  ^J  (^  (r^  Z   ■    -c  /3)  c/^  of^ 


(16) 


-r(^^f3)     =^       ^.     i^ril^j  -^     ^^  fr  l^^-^J  ^  (17) 


G(c^j-^,^}=     ^ 


^U,(e^r)UJe„^)    e    ^ 


n-  / 


(18) 
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U^(e^r)  is  th«  linear  eoabination  of  Baasel  funetiona  of  order  one 

U,  fen  rj     =    J^  (6:^r)  y;  fe„  r,)    -  J,  fe^  a/    K  {e„  rj  (19) 

and  the  eharaot eristic  y&laea    6^  are  ehosen  for  these  prjrtioular  bound- 
ar7  conditiona.     The  quantity   *^  is  the  nona  of  ^x(^j)  over  the  in- 
terrsl  r^,   r^* 

Kno\djag  now  the  radial  and  tangential  relocitiea  we  hare  froa 
Sq,  16  the  linearised  axipil  veloei^ 

^    =     ^,    -  I  -^  f-^  (ruj    dj3  (16) 

Ve  have  fonnilated  alx)ve  the  linearised  solution  as  developed  ligr 
Marble  •     It  should  be  noted  th;>t  the  Oreenis  function  as  derived  is 
independent  of  the  Banner  of  linearisation  provided*  of  coarse,    that 
the  unhor^geneous  part  f(^,^)  is  known.     If  then  ve  hare  another  value 
of  f(^«/^  ),   more  exact  then  the  above  linearized  value,   we  can  use  the 
sane  integr<0.  (iSq*  16)  to  determine  a  more  exact  radial  velocity  dis- 
tribution. 

This  inte^al,   involvin£;  an  infinite  subi  of  Beasel  functions,   is 
extremely  difficult  to  calculate.     A  solution  for  a  particular  blads 
loading  was  obtained  by  Karble  •     Tor  the  axial  flow  machine  where  the 
boundaries  exclude  the  rei^ion  near  the  axis,   i.e.,   the  small  values  of 
the  ar^inents  of  the  Bessel  function*  we  can  use  the  asymptotic  values 
of  the  Bessel  function  with  good  accuracy.     This  makes  it  possible  to 
sum  the  infinite  series  so  that  the  Ore«n*s  function  can  be  expressed 
in  closed  fom. 
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c)  Aaymptotlc  Value  of  C?reon'o  F;yictiOR 

An  asymptotic  approximation i  in  olosed  form, is  desired  for  the 
Green's  function 

where 

VJ^nr)  -   J,(e^r}  K(^nr,)'J,(e„r,)/,(€:^r]  )  (19) 

J   z         _         rj-  Wo  (^n  ra)-  n^  l/^  (en  J  _  (20) 

and  ^f,  i3  found  froa  tho  boundary  oondition  that 

O.(e^r^)-    J,(enr^)  Y/e^n]  -  J,t6nr,)K(^^rJ  ^O  (21) 

For   '♦large*   (c^r)  asyraptotic  raluos  of  tho  Dossel  functions  are 


found  in  Jahnke-ande   . 


377  , 
COS  ((:„r-  V^; 


<y,  t:^o  r) 

«.x^ 

y^7r6:„r 

Oo  (tf,r} 

v.^> 

/t  /re„r 

r,  a„  rj 

.y^ 

Jfrre„r 

K>  f^r,  H 

~^ 

(22) 


-15- 
The     coriespondlng  oharacteristio  Tolues  €-,^are  found  from  the 
solution  of     Bq.  21,  with  the  fleaael  functions  replaced  by  their  aaymp- 
totic  VBluns,   :>iA.  22. 


G 


n  =  ;r-:7f   ;   ^  =  /.  ^y  ^^  ■"-  (23) 


Uslnc  thei^e  oharacteristio  values  in  rkj.  22  we  find  the  canplete   asymi>- 
totio  value  of  the  Bessol  funotions,  caiipctlblo  with  the  boundary  con- 
ditions of  the  problem. 


J,  (^^^1      ^ 


COS 


/  nrrr        377  \ 


/  n  rr  r-  77-] 
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The  asynptotic  value  of  U  ,   U,.   and  i)^  follow  easily  and  the  esynptotio 
value  of  the  Ciroen's  function  can  be  written  as  the  Infinite  sum 

Ttxis  series  can  be  suriaed  and  expressed  in  clooexl  fonn. 

, ,        2-/3  (r-r,)'f-(<^-r,) 

^  Gfr,.;  ^.,)  ^^/Z  /„     --^"^^---^"^-r-7—  (25) 

'^^  ^  .L         2'^  rr.rJ-('<-r,) 

I  cosn  rr     ■     ■  ~  cosrr  

'/^  r,-r,  r^   -r, 
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d)  Tho  Itoration  rrocodure 

Uslne:  ^^o  as^Tiptotio  ralue  of  tha  Green's  funotion  as  dei'lrod  in 
%\q  previous  soot ion  throe  general  equfltions  can  be  written  for  the 
iteration  procedure. 

!ia  will  uae  tho  velooities  result in(3  from  a  solution  of  the  lin- 
earized problem  to  obtain  the  second  approximation  and  will  obtain 
subsequent  approximations  using  relocities  of  the  preceding  approxima* 
tiom*     Hare  we  are*   in  effeot,   assuming  that  the  rortioit/  is»   in  each 
ease*   transported  by  the  relocities  of  the  preceding  approximation* 

The  radial  and  tangential  force  components  are  prescribed  and  we 
start  with 

%=    0 

w^  =  mean  axial  velocity 
o 

If  the  velooities  w2iich  are  to  be  obtained  by  the /7th  approximation  are 
denoted  by  the  subsortpt  n  the  equations  for  tlio  /7th  approximation  vel- 
ocities are 


/3r^    eO  o*  r/^ 


/f±!^LfL  -  ^n..(n^^  '^-^]d^      (26) 


^«       =       /  /  -/"       /'^    n)   r.^^  .  .  ^    a)   ^^  ^ /^  (27) 


/3-  -oo     oi^r 


idiere 


rosh  ^#-4   -  cos  rrVl^ — 

^  />^ -I  •  ^   a\  -        I       /«<      /  ri-r,  r t.  -  r , 


(r^i-^^fi)-,    -S-^    /„ 
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f3  ^   Z 


-/     i^h'^'^^  ^^ 


^^    =     ^-o     -y  7r.S7^^^W   ^/^  (28) 

/5=  -  — 


These  three  equations,  26,  27,  28,   are  the  basis  of  the  iteration 
procefla  which  is  adapted  to  ncchanical  calcolations  in  the  next  section, 
v/e  will  wake     use  of  the  inportant  fact  th^t  0(r,2;  o^^/s)  is  the  sane 
for  each  iteration. 
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lY.  Adaptation  of  the  Iteration  Prooaas  to 
Maohanloal  Calculation  Methods 

The   integrands  of  the  expreaaiona  for  the  tangential  and  axial 
velocity  canponenta,  Uria.  26  and  28,  are  regular  over  the  regiona  con- 
sidered and  can  be  intesrated  by  any  nienerical  method  with  relatively 
little  difficulty*  Ivlechanioal  calculations  can  be  porformsd  usine  the 
equations  in  their  present  form. 

The  integration  of  the  equation  for  the  radial  velocity  Sq«  27 
beoooes  rather  involved  by  any  method,  and  for  mechanical  calcula- 
tiozis  the  form  of  the  equation  must  be  radically  changed*  An  inves- 
tigation of  the  Green's  function,  0(r,z;  o^,/s)  will  reveal  a  logar- 
ithmic singularity  at  t<=  r,  /?=  z*  To  determine  the  radial  velocity 
at  any  point  rtS  the  inteoration  must  be  carried  out  over  the  entire 
region*  The  logarithmic  singularity  of  the  integration  represents 
the  influence  at  r,z  of  an  axial  change  of  vortioity  at  r»z* 

It  will  be  advantageous  to  perfora  the  integration  in  three 
parts,  one  part  being  over  the  regiUar  z*egion,  away  from  the  singular- 
ity, and  the  other  two  parts  bcinc  over  a  ffnall  region  which  includes 
the  singularity.  Integration  over  the  latter  region  containing  the 
singular  point  consists  of  evaluating  a  regular  part  and  evaluating 
a  part  containing  a  purely  logarithmic  singularity,  ^th  thia  in 
mind  ue  write  the  radial  velocity  as  the  sum  of  three  integrals* 

U  ^  u'  ^.  u"-H  u**'  (29) 


-19- 
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a 


2-e     r-S 


f^^.^^f;   /n   i^^J.{-^f    J^dp      (3.) 


Obviously  the  integrand  for  u*  la  finite  but  we  need  to  show  that 
the  integrand  of  u**  is  regular*  It  is  easily  shown  that 
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There  is  a  slngularlt/  in  this  integrand*  however,  but  it  does  not 
ooour  near  ^^0  =   rtZ*  The  distance  from  r»z  to  the  singularity^  nieas' 
ured  radially,  is  siren  by 


-  -^(r;-?.) 

r-r-,           ± 
,     0<   r^.r,    <    2: 

^    ^('-7^.) 

f           f-r. 

From  tliese  eq.ua tions  this  singularity  is  seen  to  lie  outside  of  any 
region  whose  center  is  at  r.z.  The  integration  for  u*'  therefore  does 
not  involve  a  singularity. 

The  integrand  for  u*'*  obviously  possesses  a  purely  logarithmic 
singularity  at  ^  /9-r,a* 
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a)   Svaluatloo  of  the  Xate^^al  for  }\\ 
Uflin^  the  notation  that 


z-e      r^  °^      /J 


Z^6       r-^        Z/ft-      ^2 

}    -  J  J  ^  f  f ' 


(3^0 


ve  have 


U 


-11      -f('^.^)Gfr2,^^^)d<:^d/^ 


^35) 


The  fimction  f  (-, ^)  will  be  different  for  each  iteration  but  the 
Green's  funotion,  C(r,z;  «^,/^  ;  will  be  the  same* 


/y  /  /  /y  /,/  /  //////./ 
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////////////////// 
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Fig.  2 


SubdiTision  of  Flofw  Field 
Into  Sknall  Reotangles 
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If  wQ  dlTlde  the  region  into  amnll  reotangles  as  indicated  in 
Fig.  2  we  oan  asau^ie  with  good  aoouraoy  that  t{°<^^r3  )  oan  be  repre- 
aented  in  one  of  these  sufficiently  ssall  regions  n  by  the  second 
degree  surface 


/^  f-^^  /3j   -   c,  ^  Cr,;,  -^  ^  c^,  /s  ^  c^^^r? 


(36) 


where  four  values  of  f ( "^^  /^  )  in  the  region  n  are  necessary  to  eval- 
uate the  constants  C.  Suppose  the  foiir  values  of  f(  «y/^  )  are 
t^-^   -—  f  ,  .  liubatituting  these  values  in  Uq.  36  ire  have  foijr  equa- 
tions froQ  which  the  constants  3  oan  be  evaluated* 


•fm.        -       ^^1      "^      ^^z     ^^  m.       -^     C ^^  f3 nt      y     C/7<»   <=<m  ftnz 
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Fjur  conctr_nt3  '^^•^   — -  C  .  sjust  be  detoininod  for  each  rectangle. 
It  Kill  bo  comrenient  to  uc©  f(  ^=<jP)   ovalnQtcd  at  the  corners  of  the 
rectangle  for  this.  Uains  the  above  notation  we  can  write 

and  can  further  define 


(40 


so  that 


n 


^ff/  J  hi    "^  C  nz   '/7Z  ^  ^/>J  •'^J  '*  ^04-  In4J 


U'     =        :ZL.        C^,In.     ^Cnz    fnz    ^^^3^3     -^Cn^In.]  (U) 


Hare  the  0  ,  depend  only  on  the  function  f  ( «<^  /^  )  and  must  be  deter- 
mined for  each  iteration.  T^e  I  .  depend  only  on  the  Green's  function 
and  can  be  determined  once  and  for  ell.  It  should  be  noted  howerer 
that  the  I^  are  different  for  each  r»z}  that  ist  we  must  have  a  cer- 
tain set  of  I   to  use  in  determining  the  velocity  at  each  point. 

m 
If  we  denote  by  u*  the  velocity  at  the  center  of  the/r7th  rectangle 

and  use  the  superscript  m  to  denote  the  corresponding  I^t  then  the 

velocity  u'  at  the  center  of  the /nth  rectangle  is 


•2^• 

The  I^  are  nu^neroue  and  are  t-iuite  difficolt  to  calculate.     How 
crer  since  the  Green's  function  was  dertred  by  oonsiderlng  only  the 
fundflBziental  equ?.tloiis  ejid  the  boundaxy  oondltion  and  henoa  is  ooropleta- 


ly  Independent  of  the  presorlbed  blade  loading  or  the  blade  ahapot 
m 
'ni 


m 
these  raluea  I  .   can  be  used  for  any  axial  flow  problem  for  which 


m 
the  hub-tip  ratio  ia  the  aeme  aa  that  for  which  the  I^  are  calcu- 
lated* 

Considerable  aimplification  of  the  problem  of  calculating  the 

a 
Iq^  ia  possible  because  of  the  ayinaetry  of  the  Careen's  function*     This 

function  is  8y!]inetriCAl  about  ^-  Zt   a  fact  which  reduces  the  number  of 

m 
I  .   by  almost  one-half*     furthermore  it  possesses  another  sort  of  sym- 
metry in  /3  boSl  %  in  that  the  influence  of  the  rortioity  at /^  on  the 
Telocity  at  s  is  the  same  as  the  influence  of    the  vorticity  at  i  on  the 
velocity  at  /3  .     <nii8  condition  i^ain  reduces  the  number  of  I  "J   by  al- 
most one-half*     This  latter  "reciprocal*  relation  can  be  used  with  the 
radial  distances  o<  and  r  for  further  simplification,  but  the  relation 

f  mat 

is  slightly  different  bocouse  of  the  factor  Jyr  in  the  Green's  function* 

The  I  -  will  of  course  depend  on  the  number  and  size  of  the  rec- 
tangles into  which  the  rogion  is  dirlded*  The  rectangles  should  be 
small  in  the  ricinity  of  the  region  where  the  blade  forces  act  and 
can  be  larger  away  fron  the  blades*  Accordins  to  the  linearized  results 
given  by  Marble  it  appears  that  in  most  cases  it  is  not  necessary  to 
carry  the  integration  more  than  about  fire  blade  span  lengths  (r2-r|  ) 
upstream  and  downstream  of  the  blades t  but  this  range  may  not  be  suf- 
ficient for  the  higher  approximations  and  it  is  beliered  adTisable  to 
calculate  the  I^  for  a  somewhat  greater  range* 
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.m 


The  expression  for  u*     is  best  expressed  in  matrix  form  for  zneoh* 
anioal  oaloulations*     We  define 
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Using  this  notation  we  can  write 


/ 


r 


W 


m  ^ 


[cr,][ir.]        M 

The  four  column  niatrioes  |Cjjjlnuat  be  evaluated  for  each  iter- 
ation. The  four  square  matrices  1  Ij^  can  be  determined  once  and  for 
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all  as  soon  as  the  aubdirlaion  of  the  field  into  aciall  reotaoglea 
io  decided  upon* 
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b)   gyaluatioii  of  the  Integral  for  u * ' 

ilie  integral  u'*  is  the  z*egular  part  of  the  integration  orer 
the  small  region  containing  the  point  at  which  the  velocity  is  to  be 
calculated 


Z-^er      /-  -tS" 


z-4r     r-r 


^f^Aifr  In 


I  coi^rr^^  -COS77      ^^.^, 
cosh  rr cw  77  '       ■'    '  —      ' 


d^^d^       (31) 


We  hare  fron  I3q.  33  the  ralue  of  the  integrand  at  the  center 
point  of  the  region,  i.e*i  at  •(^^-r^zi  aa 

It  oan  be  shown  that  this  point  is  a  saddle  point  of  the  logarithmic 
function  in  the  integrand  since  it  represents  a  maximum  with  respect 
to  o<  and  a  minimum  with  respect  to  /3  • 

It  follows  then  that  if  the  region  is  sufficiently  onall  and  if 
the  function  f(  °<  %f^)   is  smooth  we  can  use  the  value  of  the  logar- 
ithnic  term  at  the  center  as  an  approxinate  niean  value  (constant)  so 
that  an  approximate  value  for  u**  is 

u"  ^    S-  In  Jsi^fl^l^h,  f  f   f^^^j  y^^^       (48) 

z-t,     r-S 
Since  this  integral  u**  represents  the  effect  of  only  one  small 

region  the  above  expression  is  sufficiently  accurate  unless  the 

change  of  vorticity  in  this  region  is  much  greater  than  anywhere  else 

in  the  flow  field. 
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^9  can  write  the  integral  ebove  for  the  /nth  reotaogle  aa 


u 


m 


-  J'"  D  (W) 


where 


,m       _      X     /^     /-  COS  IT/y^^ 

la  eTAluated  at  the  /nth  point  and 

D"^     -     I      i       'ff'-^f^j  d^d(3  (51) 

Z-fe       r'  -  f 

la  taioen  orer  the  /nth  reotangle* 

Ttie  J  are  determined  onoe  and  for  ell  aa  aoon  aa  the  aubdiriaiona 
of  the  field  are  decided  upon*  Ibe  E^  must  be  determined  for  each 
iteration* 
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c)  Svaluatlon  of  the  Integral  for  u'** 

The  interral  u"*  is  the  singiilar  part  of  the  integration  over 
the  sciall  region  oontaiziing  the  point  at  whioh  the  relooity  is  being 
oaloulated. 


2-fr     r- 


-«r 


li'e  define  a  new  function  roamentarily  as 


and  writ©  the  integral  as 


LL'" 


The  function  gi^*/^  *r)  is  regular  and  possesses  higher  order 
derivative     so  that  it  oan  be  expanded  about  the  center  of  the  rectangle 
in  a  ^claurin*s  series  in  the  two  variables  <^  •  /^  • 


-A 


(34) 


If  we  substitute  this  expansion  of  g(o<t/3)  in  txhe  integral  u" 
and  neglect  ternis  of  fourth  and  higher  orders  we  have 


u"'  -  ^(c  ^)  K,  ^  f(mr,.  K^  -  ^fM.  ^^      ''"' 
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where 
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T^e  integrals  involving  the  odd  pavers  of  (^-r)  or  (/'•s) 
Tanish  beoause  of  the  synmetry  of  the  logarithmio  teirn*     The  above 
integrals  oari  be  evaluated  by  straightforward  integration  by  parts 
to  give 
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Tte   three  constants,  K^,  yi^*   K^i  are  determined  onoe  and  for  all 
as  soon  as  the  Bubdivlsions  of  the  field  are  decided  upon* 

By  straightforward  differentiation  we  can  express  the  derira- 
tivea  of  g(  o<  ,  J  ,r)  in  terms  of  f(*<'t/^  )  and  we  can  write  the  integral 
u***  for  the  velocity  at  the  rectangle  m  as 

^here  il ,  -'I  .  K-  are  defined  above  and  where 
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d)  The  OompXete  Integrid  for  u 

vve  have  oraluated  the  integral    for  the  radial  Telocity  in  three 
parts  as  given  by  Eqs*  4^t  ki*  ^^  39*     Oooibining  these  results  we 
oan  write  a  oomplete  expression  for  the  radial  velooity  at  the  center 
r,«  of  the  rectangular  subdivision  m. 


Z/  '^  "^  /77  /T7      / 

The  constants  C»  D  and  I^  are  eraluated  from  the  unhaoogeneous 
pert  f(«<t/9)  of  the  differential  equation  for  the  radial  velocity 
(lilqs*  11  and  27)  and  are  to  be  determined  in  the  manner  described 
in  the  preceding  seotioast  These  constants  must  be  detenoined  for 
each  iteration. 

The  constants  I«  J  and  K  depend  only  on  the  Oreen's  functiont 
(i.e«»  the  boundary  conditions),  the  hub- tip  ratio t  and  the  size  of  the 
subdivisions  and  can  therefore  be  evaluated  once  and  for  all*  I^ey  are 
used  in  the  s&iie  manner  in  each  iteration  step* 
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T.   DI3CUSSI0H 

A  Bethod  of  soItId^  the  h7drod7namic  •qoations  for  the  InooB- 
preselble  flov  of  an  inrlseid  fluid  through  an  axial  flow  turboBnaohine 
has  l>e8n  dereloped  as  an  iteration  procedure.     The  next  logleal  step 
of  this  approach  to  the  prohlea  is  the  evaluation  of  the  invariant 
terns  of  the  iteration  expression  for  the  radial  velocitijr.     Once  this 
is  done  for  the  several  likely  hub- tip  ratios  it  will  he  relatively 
easy  to  ealculite,  by  meehaniCAl  means,    solutions  for  any  axial  flov 
machine  v;ith  any  prescribed  blade  loading. 

Zn  detexmining  the  degree  of  acouracy  th^t  vill  be  required  the 
restrictions  (incorapressibility,   etc*}  which  were  iaposed  to  siiaplify 
the  nathffinatical  problem  cnst  be  considered  and  evaluated.     The  solu- 
tion should  be  consistent  vith  these  restrictions  and  of  sufficient 
accuracy  to  indicate  the  r-roprr  trends  of  the  variables.     It  is  not 
certr.in  thr.t   the  line&rlxed  solution  lacets  these  requirements  in  all 
oases.     The  iteration  procedure  can  be  used  to   detenaine  the  accuracy 
of  the  linearized  solution  and  if  necessary  to  obtain  solutions  of 
greater  accuracy. 
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VI  ATFMDIJC 

a)  A  Brief  Conp^rlaon  of  the  Axial  riov  and  Hlxod  now  FroUeai 

The  extension  of  the  nore  recent  axial  flow  tolutiona  to  apply 
to  mixed  flow  machines  vould  be  extrcnely  difficult  and  vould  rec^uir* 
consideration  of  sarer&l  points,  not  significaoit  in  axial  flow,   vdiioh 
are  of  ataost  importance  in  nixed  flow.     Three  essential  differences 
between  the  two  problems  are  x>ointed  out  here  in  order  to  indicate 
which  of  the  assunptions  used  in  the  axial  flow  analysis  would  not 
be  applicable  to  a  mixed  flow  analysis. 

(1)  The  nixed  flow  machine  contains  continuous  Tanes  as  con- 
trasted with  the  rotor  and  state r  blade  rows  in  an  axial 
flow  Q^jchine.     The  ranes  therefore  cannot  be  "twisted** 
without  introducing  excessiTe  tilt  away  from  radiel  so 
that  the  oentrifu^sl  forces  cause  large  bending  moments 
in  the  vane  and  prevent  operation  of  the  machine  ^t  ex- 
trenely  hi^  speeds.     A  geonetric  relation  oast  hold  be- 
tween the  relative  velocity,   the  vane  foroes  on  the 
fluid  Mid  the  shape  of  the  vanes  throughout  the  region 
where  the  vanes  are  present,  ^^lereas  for  blade  rows  as 
in  axial  flow  this  relation  holds  only  in  the  region  of 
the  narrow  bl£de  snd.   In  fact,   may  be  concentrated  in  a 
"lifting  line**  for  a  good  aiTriroxiiaation, 

(2)  The  lajTger  radial  velocities  which  naturally  occur  in 
mixed  flow  prohibit  the  simplifying  assumption  used 
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for  axial  flow  thnt  the  radial  Telocity  Ib  rezy  small 
oonipj*red  to   the  axi-^l  velocity. 
(3)       The  boundary  conditions  of  the  aathea^tlcal  problem 

of  the  cilxed  flov  oomj>res8or  are  sxuch  more  coBzpllcated 
than  those  for  axial  flow  end  gre&tly  increase  the  dif- 
ficulty of  obtRlning  a  solution.     The  Tarlables  are  not 
separable  in  this  case.      It  should  be  noted  too  that 
the  boundaries  will  probf.bly  be  very  different  for  each 
mixed  flow  problem  whereas  they  are  always  essentially 
the  sasie  for  axial  flow. 
Several  types  of  Tones  are  x>osalble  but  two  special  vane  shapes 
are  likely  to  be  of  interest. 

The  first  of  the^e  might  be  called  "radial  Tanes".     These  raaes 
are  generated  by  redial  lines  through  the  axis  of  rotation.     Here  the 
radial  force  is  zero.     Badial  Tpnes  are  necessary  In  a  high  speed 
ffl&ohine  because  of  the  high  centrlfugpJ.  forces.     If  radial  velocities 
and  (or)  pitch  angles  are  large  as  in  mixed  flov  then  the  angles  be- 
tween the  vanes  -nd  the  hub  or  shroud  will  be  acute,    thus  increasing 
bound&iy  layer  effects. 

The  second  type  might  be  called  "normal  vanes'*.     These  vanes  are 
generated  l3y  lines  through  the  axis  but  tilted  in  a  meridional  plane 
so  AS   to  be  noxQal  to   the  meridional  trace  of  the  streanllne  at  all 
points.     Here  the  vane  force  has  two  cooponents,   one  acting  alon^  the 
meridional  trace  of  the  streamsurface  and  one  acting  t?^ngentially. 
TMs  vane  is  not  structurally  adequate  for  extremely  high  speed  rotors 
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but  doea  hpve  the  good  tealMre  th^t  rrio^les  between  the  Tane  L^nd  the 
hub  or  shroud.  :ro  ri^^t  aneiles,    a  f^^t  w.lch  nay  tainlmli©  boundary 
leyer  effects.     It  i3  interostlng  to  note   that  one  force  component 
Accelerptes  the  fluid  In  its  path  in  a  meridional  plnne  and  the  other 
cooponent  accelerates  it  t-i^^entially  bo  th&t  it  appears  th?.t  no 
"vnated"  forces  pTe  present*     For  the  axial  flow  machine  radial  ranes 
and  nornal  TfJies  are  about  the  aa&e. 

The  selection  of  either  radial  or  normal  T<'«nea  vill  lead  to 
great  sioplificaticn  of  the  mathen-tical  problem. 
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J?)   30^3  Oomocaonceg  of  the  Hydrod:T"n.-^mi o  ^(^uations  ag  Applied  tp 

Im  2>aler  equations  of  notion  are  written  in  vector  form  and  in 
c.^lindricrol  oon^onent  fom  i'or  the  isontro -ic  flow  of  an  inviscld 
rind  incoiaprassiole  fluid  acted  upon  by  non-conserrr.tive  body  forces. 
Tha  flow  is  symmetrical  about  the  axis. 


i-  ^  ^  iA,  ^^'^  -         F^  /gj 


Or  £)z-  r   c^^ 

BQr  use  of  ?  vector  identity  the^e  equations  n?»y  be  written  in 
,|^  terms  of  tho  vortlcity 

where,  since  tho  flov  is  isentropic  and  adiabatio, 


rrf;  =  7^<^p  ^  ± 


I/- 1' 


is  the  gradient  of  the  totc-0.  nnergy  of  the  fluid.  In  component  form 
■Sq.  (3)  is 
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c^eToral  interesting  relations  between  velocity,   TortieiV*   total 
energy,   rane  ahapOt   etc,   oan  be  derived  from  these  equations  if  vt 
conaider  the  flov  between  t\^  concentric  bounding  •urfpjces  acted  npoA 
by  Trnes  or  bl.^dea  rot-tin^  with  an  an^^tr  roloci^  ^.     Ve  will 
833une  pji  infinite  nunber  of  rpjies   (or  blades)   so  thet  the  forces  of 
the  Yrjie  on  the  fluid  raay  be  rerresonted  by  body  forces  synmetrioal 
about  the  axis.     Turthemore  since  the  fluid  is  inviscid  these  forces 
are  nomal   to   the  vane  and  hence  to  the  rel retire  velocity  so  that 

l^f--  ooT-J  ■  F    ^  o  (6) 

If  've  t*>:3  the  scsl-r  product  of  the  velocity  and  both  sides 
of  Eq.   3  these  reralts 

o 

or 


^.(^]  =  ^-^^ 


(•) 


yrom  the  second  equation  of  notion 


Prora  Bqs.   6  and  7  ve  see  that 


—.  ( -T^j      =    ^  -^^ryj       ,   vhere  rotating  vanes  act 

-      O  J  ^Aiere  no  rotating  vanes  act 
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If  ve  oonslder  the  ca^e  where  the  fluid  has  constant  energy 
far  upttreaa  and  is  fre^  cf  Torticity  far  upatrean  and  make  the  stlp* 
ulatlon  thp-t  the  total  energy  of  the  fluid  ia  changed  onl;'  tj  the 
action  of  the  raovinc  vanes   (I.e.,   ad.ia"batlc  flow}  thcto  3q.   8  can  be 
integrated  along  each  streoaeiirface  with  tr.o  result  that 

for  thit  case  of  uniform  ener^  and  no  Yorticity  fer  upatrean,       ^ 
and    cDrt^  are  constant  far  upstream  so  thati 


pr /-^  y   =z  O  ;  upetream  of  the  Tanes 

—  ^((^r'^)  ;   in  regions  of  and  dovmetreaa  of  the 

ranee  if  no  forces  such  as  stationary 
Tana  forces  occur 


(9) 


If  any  fercos,    such  as  those  resuiltini:  from  stationary  Tanes  or 
stators,  hare  tnn^entir.l  eooponents  then  we  can  only  say  that   0T\r'/  ~  ^ 
alonf  the  stroaasurface  since  in  this  case     oO  'i-fryj  =  o  ,  b-ui    ^  (n^j  ^  o 

Oombining  ?3qa.   Z  and  9  we  hare 

Fa  _Q      ^     s^  (cDr\r)    '  F 
and  if  ws  Htdtiply  l)oth  sides  by  the  Tortioity 

V6  see   that 

(10) 


f2    '  f^fuir^}  -  F J 


-    a 
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TTslrif:  only  the  doflnition  of  ▼ortlGity  it  is  easily  ahovn  that 

n.   •    vfcori^)    ^  o  (11) 

Qnd  hence 

f2    •    i^-{/j     =  O  (12) 

and  further  froa  7'^,  10 

22  ■  F       ^o  '  <") 

Using  the  definition  of  vorticity  liid  3q«   9  it  is  seen  that 

cD?  .  C2     -    ^  f^)  (14) 

Combining  ^q,  14  and  Tlq,  3 

?or  this  special  case  "re  can  draw  the  following  interesting  con- 
clutions: 

Upstrean  of  the  Vines i 

1.  The  Torticitv  Is  zoro, 

2,  The  totpJ.  en.^rgy  of  the  fluid  is  uniform. 
In  the  Region  Miere  Yane  Forces  Actt 

1,  The  vorticity  rector  is  tangent  to  the  vane  surface.    (Sqa,  5 
and  13). 

2,  The  vorticity  vector  is  tangent  to  a  surfrce  of  constant 
total  energy.      (Sq,  12) 
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0.  riie  uound  vorticity  is  tingent  to   the  lir.j  of  intersection 
of  £  v&ne  snd  a  concentric   aurf^jcs  of  Cjaataiit  tot.l  eaargj' 

■J'Oynstreii:;  of  tl;c  Vr.r-Oij; 

1,  L':.:  vort-lc'lty  v-ctor  is   tr'nc-^r^t   to   tlit-   rcl:,tive  VGlocity 
vector  (reltlve  to   t!:*.*  rot  tin<':  re  tor}, 

;. .  ?'.6  vorticity  voctor  1?   tor-gent  to  ?  z\xrt?.ce  of  conatent 
totpl  enerr^v',    this  surface  being  the  gtreajr-siirft-co. 
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